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We compare recent, seemingly different, approaches to TMD-factorization (due to Echevarria, 
Idilbi, and Scimemi and to Collins), and show that they are the same, apart from an apparent 
difference in their definition of the MS renormalization scheme. 



I. INTRODUCTION 

In a number of processes, notably the Drell-Yan pro- 
cess at small transverse momentum, it is important to 
use transverse- momentum-dependent (TMD) parton dis- 
tribution functions (pdfs) and an associated TMD factor- 
ization theorem of QCD. However, there is a variety of 
different formulations in the literature, with seemingly 
very different methodologies and different definitions of 
the TMD pdfs. 

Some of the differences even amount to apparent in- 
compatibilities. For example, in the formalism developed 
by Collins, Soper, and Sterman (CSS) [HE], ™ D P dfs 
depend not only on the essential kinematic parameters 
x and kr, and on a renormalization scale /i, but also on 
an extra auxiliary scale C s . An essential role is played 
by the Collins-Soper (CS) evolution equation [1 of the 
pdfs with respect to £. In contrast, in the recent work 
by Echevarria, Idilbi, and Scimemi (EIS) [3] in the con- 
text of soft-collinear effective theory (SCET), there ap- 
pears to be no corresponding parameter. Indeed, in a 
phenomenological application of this work by Echevar- 
ria, Idilbi, Schafer, and Scimemi [4], there even appears 
the statement that they "do not have to solve Collins- 
Soper equation". Furthermore, in the recent improve- 
ment of the CSS method by one of us (JCC) in Ref. [5], 
the presence of the £ parameter appears to be intimately 
associated with the use of non-light-like Wilson lines in 
JCC's gauge- invariant definition of the TMD pdfs. But 
all the Wilson lines in EIS's definition are light-like (with 
the aid of their ^-regulator in a limit that the regulator 
is removed). Finally, the definitions in the two meth- 
ods appear very different: In the JCC method the finite 
TMD pdf is [5] Eq. (13.106)] a renormalized convolution 
product of a bare TMD pdf with a square root of three 
bare soft factors, two of them with the characteristic non- 
light-like Wilson lines. But in the EIS definition [3] Eq. 
(2.13)], the square root involves only one bare soft factor. 
JCC and EIS do agree on the presence of the square root. 



In this paper, we show that nevertheless the EIS and 
JCC definitions give exactly the same TMD pdfs (and 
hence the same hard scattering factors in TMD factor- 
ization properties). We show how the EIS definitions do 
in fact have a parameter £ that EIS have identified with 
Q 2 . Thus their TMD pdfs obey a CS equation. Although 
the application of the EIS work in j4] does not explicitly 
use the CS equation, it actually uses equivalent results, 
the most notable of which are the results of Korchemsky 
and Radyushkin [6] on the Q-dependence of the anoma- 
lous dimension of a cusp anomalous dimension. We show 
how the JCC definition of the TMD pdfs can be cast into 
exactly the form of the EIS definition, in which no non- 
light-like Wilson lines appear and in which there is only a 
single bare soft factor instead of three. It should be noted 
that, in both cases, the bare factors are separately diver- 
gent because of the presence of light-like Wilson lines, 
which must be regulated, and a limit taken with the reg- 
ulator removed. The EIS and JCC differ in the preferred 
method to regulate light-like Wilson lines. The EIS form 
of definition requires coordination of the rates at which 
are removed the regulators of different light-like Wilson 
lines. 

EIS do say in [4 that "one can argue that the two 
[methods] are equivalent", but they do not provide any 
details of the argument Q 

We expect that the methods developed in this paper 
can be applied to analyze and/or relate many other meth- 
ods for formulating TMD factorization. A list of some 
relevant papers beyond those cited above is [8HT4]. 

The reader needs to be very aware that the JCC 
method we use here differs in very important details 
from the earlier CS formalism on which it is based, and 
also that there was a change in candidate definitions of 
collinear factors in the course of the derivations in [5]. 
The definitions that we use are of the form given in [5j 
Sec. 10.11] for the Sudakov form factor, and in [5, Sec. 
13.15.3] for TMD pdfs (with a modification to use past- 
pointing Wilson lines for the Drell-Yan process, as sum- 
marized at the beginning of [5[ Sec. 14.5]). 
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1 See also [7], which appeared as we were completing this paper. 
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The paper is structured as follows: In Sect. [II] we in- 
troduce notation and enumerate the common features of 
the EIS and JCC methods. In Sect. Ill we prove that 
the EIS and JCC methods both give the same result for 
the cross section. We prove in Sect. [IV| th at the JCC 
and EIS TMD PDFs are equal. In SecT|V| we identify 
the parameter in the EIS TMD PDFs that corresponds 
to £ in the JCC method. We summarize our results in 
Sect. Enl 



II. COMMON OVERALL STRUCTURE 

We work with the Drell-Yan process: production of 
a high-mass lepton pair of 4-momentum q from a sin- 
gle virtual electroweak boson in a collision of hadrons of 
momenta pa and pb- The invariant mass of the lepton 
pair is Q = ^q^q^, its rapidity is y, and its transverse 
momentum is qr- 

When qT <C Q, TMD factorization as formulated by 
both EIS and JCC is that the cross section is proportional 
to a Fourier transform in transverse coordinate space of 
TMD pdfs: 



da 
¥q 



J d 2 b T e 



ibT-q T 



fi/A&A, b T ) fi/B(xB, b T ) Hij(Q) + p.s.c. (1) 



Here, the sum is over the relevant parton flavors, fi/A 
and fj/B are Fourier transforms into transverse coordi- 
nate space of TMD PDFs, Hij is a hard scattering, and 
"p.s.c." denotes "power-suppressed corrections". The 
corrections to the formula are either suppressed by a 
power of qT IQ or by a power of M/Q, where M denotes a 
typical scale for light hadronic masses; those corrections 
suppressed only by a power of qr/Q can be handled by 
normal factorization methods with integrated pdfs, as by 
the Y term in the CSS formalism. 

To simplify the notation, we write the above structure 



a = ABE, 



(2) 



with all the kinematic variables, flavor indices, and cor- 
rections left implicit. The intuitive interpretation is that 
the factors A and B contain dependence on momenta 
collinear to each of the corresponding beam particles. 

As for the kinematics, we have xa%bs = Q 2 and the 
rapidity of the Drell-Yan pair in the center-of-mass frame 
of the beams is y = \ In In light-front coordinates, 
the beam momenta are 



Pa 



Pb 



2p~ L 



,P b ,0t), (3) 



with p\ and p^ large, so that s ~ 2p\p^. The plus 
momentum of the parton from A is k\ = xap\ , and the 
minus momentum of the parton from B is kg = xbPb- 



The above statement of TMD factorization is used by 
both the JCC and EIS methods. However, in we have 
indicated only the dependence on kinematic parameters 
of the factors fi/A fj/B, and H^. In addition all three 
factors depend on auxiliary parameters. It is common to 
both methods that all the factors depend on the renor- 
malization scale ja. There is one further parameter for 
the TMD pdfs in both methods. 

In the JCC approach, each collinear factor also de- 
pends on a parameter. This is constructed from an arbi- 
trarily chosen rapidity parameter y ni which has the intu- 
itive purpose of separating right and left moving quanta. 
For A, the £ parameter is 



( A = 2(fc+) V 



2j/„ 



2(x AP +) 2 e- 2 y", (4) 



and for B it is 



CB=2{k- B fe 2 y-=2{x BV - B ?e 2y -. (5) 

These two quantities obey (a(b = Q 4 - In both the EIS 
and JCC approaches, neither the hard scattering H nor 
the cross section has dependence on the ( variables. 

In the case of the EIS method, the authors' statements 
and notation suggest that the TMD PDFs depend on 
H but not on some other parameter like £ — see the 
paragraph above [3, Eq. (1.6)], and see [3[ Eq. (2.15)]. 
In fact, the EIS PDFs do depend on Q as a separate 
parameter distinct from \i. This can be seen from [3j 
Eq. (2.17)], where there is Q and \i dependence for the 
coefficient C n that relates the TMD PDF to the ordinary 
integrated PDF. It can also be seen from the fact that 
the anomalous dimension of their TMD PDF depends on 
both Q and \i — see [3j Eq. (3.21)]. In later sections 
of this paper, we will show how the Q 2 parameter of 
the TMD pdfs with the EIS definition should in fact be 
identified with the £ parameter of those with the JCC 
definition. 

Common to both schemes is the property that the 
product of the two TMD functions is related to a combi- 
nation of unsubtracted TMD functions and a soft func- 
tion: 

AB — ^°J^° x UV renormalization. (6) 
So 

Note that all the authors agree on applying UV renormal- 
ization in the MS scheme to remove divergences where 
transverse momenta go to infinity Such divergences 
are all regulated by dimensional regularization. 

Each of the unsubtracted functions is defined with the 
aid of light-like Wilson lines, taken as a limit ^\ The Aq 



2 But see Sec. |Vl| for a possible difference in the definitions of the 
MS scheme. 

3 Note that these Wilson lines must be joined by transverse links 
at infinity to give exactly gauge- invariant quantities |15l 116) . 
Only with these transverse links can correct results be obtained 
in all gauges. Notably, the transverse links are essential to get 
the correct Sivers asymmetry in light-cone gauge ^] See 
Refs. |17l I18| for corresponding statements in SCET. 
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and B are the "unsubtracted" TMD pdfs / unsub in [5j 
Eq. (13.106)], or the f n functions in [3 Eq. (3.21)]. The 
So is a light-cone regulated soft function, defined accord- 
ing to [5, Eq. (13.106)] and [3, Eq. (2.7)]. The division by 
So cancels double counting of soft gluon contributions be- 
tween the Ao and Bo factors and also cancels the rapidity 
divergences that would otherwise occur because light-like 
Wilson lines are used in Ao and Bo- Since each of the 
factors in ([6| separately contain light cone divergences, a 
meaning can only be given by applying regulators to the 
Wilson lines, and then taking the limit of removing the 
regulators in a manner consistent with factorization. 

As explained in [5 , Sec. 10.11.1], the limit of light-like 
Wilson lines and the removal of dimensional regulariza- 
tion do not commute. We will always define that the 
limit light-like Wilson lines is to be applied first. Then 
MS renormalization is applied and dimensional regular- 
ization removed. EIS's work appears to use this order as 
well. 

At this point, the only relevant difference between the 
EIS and JCC treatments is the choice of method for reg- 
ulating the divergences that come from using light-like 
Wilson lines. The JCC method uses non-light-like Wil- 
son lines, and takes the limit that they becomes light-like: 



Ajcc^jcc 



lim 

VI — S- + 0. 



AjCC,o(^2)ffjCC,oQ/l) 
^00,0(2/1,2/2) 

x UV renormalization. (7) 



(This result follows directly from [5, Eqs. (13.106) and 
(14.31)].) The Wilson lines appear in Feynman diagram 
calculations as eikonal propagators of the form, 



(8) 



i0 7 



-e~ 2 yik+ + k~ -M0' 



together with the appropriate vertex factors. Here, y\ 
and y2 are the finite Wilson line rapidities, and the signs 
are those appropriate for when the gluon momenta are 
oriented as in the one-loop graph for the Drell-Yan pro- 
cess in Fig. [T] Note that, in the JCC method the reg- 
ulated factors are all matrix elements of exactly gauge- 
invariant operators. 

The EIS approach [3] uses what they call a ^-regulator 
instead of ([sj) , where the Wilson line propagators for reg- 
ulated light-like lines have the following forms: 



(9) 



-k+ + i5+ k~+i5~ 
Their equivalent of Q is 

a o r ^EIS,0(^ + )^EIS,0(^") 

^eis^eis hm — — — 

5+,5-^0 5eIS,o(£",^ + ) 

x UV renormalization. (10) 

That is, instead of using Eq. (J8|, they replace light- 
like propagators with ones with fixed non-zero auxiliary 
terms inserted into the denominators. 



FIG. 1. One-loop graph contributing to the Drell-Yan process. 



The right-hand-sides Eqs. Q and (10) differ only in 
how the light-like Wilson lines are regulated, and there- 
fore the limit as the regulator is removed must be the 
same in each case. The same UV renormalization pre- 
scription, MS, is used. Hence the product of the two 
TMD pdfs on the left-hand-sides is the same, i.e., AB is 
the same in the EIS and JCC schemes. The remaining 
issue is to determine whether the individual TMD pdfs, 
A and 5, are the same in the two schemes. 

Our proof that this is indeed the case will use factor- 
ization properties for the unsubtracted factors, and one 
immediate consequence of these factorization properties 
is that the limits in Eqs. ([7]) and (10) do exist. 

One complication is that EIS also use a correspond- 
ing regulator for the quark lines in their calculations, for 
example, [3, Sec. 3]. This provides a regulator for the 
collinear divergences associated with quark-gluon inter- 
actions in massless QCD. The regulators for the quark 
lines and the Wilson lines have to correspond: Their reg- 
ulator parameters of the quark lines are denoted A + and 
A - , and they set A + = S + Pg and A - = S~p\ (in our 
notation). This coherence is needed so that the combi- 
nation AqBo/So is a valid approximation to the original 
graphs for the cross section. 

However, for defining the parton densities, and for us- 
ing ( [To] ) , one can ignore this issue, provided that we also 
use dimensional regulation. Since dimensional regulariza- 
tion regulates the collinear and soft divergences, we can 
take the limits A ± — > first. Then the limits S ± —> 
can be taken separately to define the parton densities. 
These results are compatible with the factorization re- 
sults obtained below. However, as we will demonstrate, 
the limits on £ + and on S~ need to be coordinated. 



III. FACTORIZATION PRESERVES 
FINITENESS AND REGULATOR 
INDEPENDENCE OF lim (A B /S ) 



To see explicitly that the limits in ^ and (10) exist 
and are the same, we apply factorization by the JCC 
method to express the regulated Aq, Bq and So in terms 
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of products of the finite collinear factors in the JCC 
scheme. Factorization is applied separately for Aq etc 
regulated by each of the JCC and EIS methods. (EIS 
would undoubtedly apply their factorization method in- 
stead.) 

In all of this paper, we will take for granted that fac- 
torization holds for the cross section and for the other 
situations that we consider. Our aim is to use factoriza- 
tion to relate different formulations. 

There are a couple of complications that affect the or- 
ganization of our proofs: 

• In constructing the TMD pdfs there are limiting 
operations: to remove the regulator (s) on light-like 
Wilson lines, and to remove dimensional regular iza- 
tion. These limits do not commute — see, e.g., [5j 
Sec. 10.8]. The definitions are first to remove the 
regulator of light-like Wilson lines, e.g., in Eqs. ^ 
and (10), then to apply UV renormalization, and 
finally to remove dimensional regularization. 

• Corresponding to the different steps in taking lim- 
its, there are actually three different objects that 
implement some kind of TMD pdf: 

1. The unsubtracted quantities Aq, etc, as on the 
right-hand side of Eqs. ^ and (10); to make 
these finite a regulator must be used on their 
light-like Wilson lines. Were all divergences 
absent, these would give the natural defini- 
tions of TMD pdfs. 

2. The unsubtracted TMD pdfs are combined 
with square roots of unsubtracted soft factors, 
in definitions that we present below, in which 
the limit is taken that the regulator (s) of the 
light-like Wilson lines is removed. (The JCC 
definition continues to have certain other non- 
light-like Wilson lines.) These we will call "un- 
renormalized" TMD pdfs, notated A unren , etc, 
and we will prove these are the same in the 
JCC and EIS definitions. 

3. Finally UV renormalization is applied and 
then dimensional regularization removed. The 
resulting TMD pdfs are those on the left-hand 
side of Eqs. ^ and (10), they are finite in 
QCD, and they are the TMD pdfs that are di- 
rectly used in phenomenological applications. 

There are several kinds of divergence that appear at 
various stages of using factorization. 

First are the UV divergences in QCD itself; these are 
removed by renormalization as part of the definition of 
QCD. 

There are also UV divergences in the definitions of pdfs 
(and other related quantities). Although these diver- 
gences are also removed by essentially conventional renor- 
malization, their status is rather different from those of 
QCD itself. The definitions of pdfs arise from the appli- 
cation of approximations appropriate to a collinear region 



of momentum, and the UV divergences arise from a nat- 
ural extrapolation of the approximations to all momenta. 
We will therefore characterize these divergences as "in- 
duced divergences". They cancel between the factors in 
a factorization property (as do the corresponding renor- 
malization factors). Induced UV divergences also arise 
for the same reasons in the ordinary short-distance op- 
erator product expansion (OPE). It is possible to simply 
cut-off induced divergences. But a cleaner formalism re- 
sults if one allows the divergences to occur and defines 
finite operator matrix elements by renormalization. 

Then there are the IR and collinear divergences that 
appear in Feynman-graph calculations with massless 
quarks and gluons. These would be absent if all the 
fields have non-zero mass, as can be true in a model the- 
ory. They are presumably cutoff by (non-perturbative) 
confinement in QCD. In perturbatively calculable short- 
distance quantities, like hard-scattering coefficients, it is 
useful to neglect masses to simplify calculations. But the 
IR and collinear divergences must cancel in such quanti- 
ties, precisely because they are short-distance quantities. 

Finally, there are rapidity divergences. These arise 
when TMD pdfs are defined in terms of operator matrix 
elements with exactly light-like Wilson lines. In a gauge 
theory such Wilson lines are a natural consequence of the 
expansions in small variables for collinear and soft mo- 
menta in setting up factorization; for this reason, they 
occur in both the EIS and JCC methods. Having sepa- 
rately defined operator matrix elements is important for 
analyzing asymptotic behavior in the form of factoriza- 
tion properties, and for identifying the operator struc- 
ture of universal non-perturbative objects like pdfs and 
fragmentation functions. These matrix elements can be 
studied with non-perturbative methods, and are of in- 
trinsic interest as objects related to hadron structure. 
Rapidity divergences are another example of induced di- 
vergences, i.e., divergences introduced by the approxi- 
mations used to factorize the cross section. To define 
finite pdfs, it is possible to cut off the rapidity diver- 
gences in some fashion (as in the original CS formalism), 
but it is also possible to remove them by a kind of gen- 
eralization of the renormalization that is applied to UV 
divergences. Both the EIS and JCC formalisms use such 
a generalized renormalization, implemented by multiply- 
ing unsubtracted TMD pdfs by suitable square roots of 
unsubtracted soft factors (or by some equivalent of this 
procedure). 

Induced divergences, no matter whether they are UV 
divergences or rapidity divergences, are unphysical in the 
sense that they do not correspond to regions of momenta 
that give actual divergences in unapproximated graphs 
for scattering. Nevertheless they are symptoms of im- 
portant physical phenomena. They arise when in some 
region of momentum one expands to low order in pow- 
ers of a small variable relative to a much large variable. 
Such expansions lead to the definition of quantities like 
pdfs. The natural extrapolation of the expansion to all 
values of the expansion parameter gives the induced di- 
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vergences. The range of momenta over which the expan- 
sion gives a valid approximation to the original graph 
depends on the energy and other kinematic parameters 
of the process under consideration. 

Thus while induced divergences, including rapidity di- 
vergences, are in one sense unphysical, they are also phys- 
ical in another sense that they are related (by extrapola- 
tion) to important, and even dominant, physical effects. 

Associated with the generalized renormalization of in- 
duced divergences to construct finite pdfs, etc, are typi- 
cally auxiliary scales like the renormalization scale \i of 
MS renormalization. Such a scale can be intuitively char- 
acterized as the range over which the underlying approx- 
imation is used. Since the range over which the approxi- 
mation is useful is process-dependent, evolution (e.g., by 
the renormalization group) with respect to the auxiliary 
parameter (s) quantifies how the actual cross section is 
related to particular ranges of momenta inside the am- 
plitudes for the process. 

The generalized renormalization of rapidity diver- 
gences in the JCC scheme leads to the auxiliary parame- 
ter that we have already mentioned, with the associated 
CS evolution equation. As we will see, the EIS method 
also uses corresponding properties. 



A. Ao, Bo and So in the JCC method 

In Ajcc,o? the hadron and the Wilson line of rapidity 
2/2 are widely separated in rapidity, so they factor into 
two collinear pieces: 

AjCcAVl) = ^JCC, unren. (2(^) 2 e~ 2 ^) 

Cjcc, unren, W (y n ~ 2/2) + error. (11) 

The factors are illustrated in Fig. [2] As usual, the leading 
regions have collinear parts, a soft part, and a hard part. 
The collinear parts are associated with the hadron and 



with the Wilson line of rapidity y^ and correspond to 
the factors on the right. No separate factor is needed for 
the soft part, according to the methods of both JCC and 
EIS. As usual in the JCC scheme, the collinear factors 
depend on a rapidity parameter y n: which we choose to 
be the same as in Q, although our argument can be 
generalized to deal with a different value. Notice that we 
indicate the dependence on the auxiliary parameters and 
the Wilson line rapidities and not the other kinematic 
parameters (e.g., x and br-) 

As for the hard part, it is in fact to be replaced by 
its lowest order value, which is unity. The reason is that 
at this stage, we keep dimensional regularization as we 
take ?/2 — » — oo. Now the hard scale is set by the to- 
tal energy of the process, which we can parameterize by 
k\e~ V2 . With dimensional regularization in 4 — 2e di- 
mensions, higher orders in the hard scattering are sup- 
pressed by a power l/{k\e~ V2 ) constantXe ? and therefore 
vanish when y<i — >• — oo at fixed e. (If we removed dimen- 
sional regularization and renormalized everything before 
taking yi — >• — oo, we would get differently organized re- 
sults; this is the non-commutativity of limits mentioned 
earlier.) Since we have not yet removed dimensional reg- 
ularization, we do not yet need to apply UV renormal- 
ization. Consequently, we are working with the unrenor- 
malized collinear factors, as indicated by the notation on 
the right-hand-side of Eq. (11). 

Apart from the issue of UV renormalization, the 
collinear factor Ajcc, unren. is the same as in the JCC 
version of Drell-Yan factorization. That is, it is the 
TMD pdf associated with hadron A, defined as in [SJ 
Sec. 13.15.3], but with past-pointing Wilson lines, as ap- 
propriate for the Drell-Yan process that we work with 
here. It includes factors involving the square roots of 
soft factors, and in two of these factors, there is a Wilson 
line with a finite rapidity, y n . The quantity Cjcc,w is a 
corresponding collinear factor associated with the Wilson 
line of rapidity y^ in the left-hand-side. 



We recall the definition (0 Eq. (13.106)]) of A JC c' 



^-JCC, unren (CO = lim Ajcc,o(^ 2 ) 



Ajcc,o(-°o) 



Sjcc,o(2/i - Vn) 



Sjcc,o(2/i - 2/ 2 ) Sjcc,o(2/n - 2/2) 
' £jcc,o(+oo - y n ) 

SjCC,0(+°O - (-OO)) SjccAVn - (-00)) ' 



(12) 



where the meaning of the second line is given by the limits in the previous line. The corresponding definition of 

Cjcc,w is 



Cjcc, 



2,w(y n ~ 2/2) = lim 5jcc,o(2/i, 2/2) 



SjecAVn - 2/2) 



£jcc,o(2/i - 2/2) £jcc,o(2/! - Vn) 



Sjcc,o(+°o - 2/2) 



#JCC,o(2/n - ("OO)) 



#jcc,o(+oo - (-00)) 5jcc,o(+oo - y n ) ' 



(13) 
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where the role of y[ and y' 2 is exchanged compared with the definition of Ajqc- 

Invariance under boosts in the z direction shows that the dependence of Ajqc on k\ and y n is only via the 
combination (k\) 2 e~ 2yn . Similarly, the dependence of Cjcc,w on y n and y<2 is only via y n — 2/2- The derivation of 
(11) is simply a version of the proofs of factorization for DIS and Drell-Yan cross sections that is given in Ref. [5]. 

We now apply the CS equation to obtain the ^-dependence of Ajcc (2(k\) 2 e~ 2yn ) and of Cjcc,w(yn — 2/2) at fixed 
values of their arguments. The CS equations are 



dy n 

unren, 

dy n 



-AjcC, unren(2(^) 2 e- 2 ^)K u 
Cjcc, unren, unren ? 



(14) 
(15) 



where K is the CS kernel which depends on br and /i, but not the other variables. These equations have the immediate 
implication that the product of A and C in (11 ) is independent of y n . The tilde in K merely indicates that this is the 



object in transverse position space, not in transverse momentum space. 

Let M be some chosen fixed reference mass scale. We use (14) to express Ajqc (2(k\) 2 e~ 2yn 
with 2(k\) 2 e~ 2yn replaced by M 2 : 



iJCC, unren 



(2(4) 2 



-2y n 



unren 

(M 2 ) exp 



In 



M Vn u 



in terms of its value 



(16) 



Similarly, 



Cjcc, unren, W (yn ~ 2/2) = CjGC, unren, W 

(O)exp \(y n -y 2 )K, 



(17) 



Applying these results to Eq. (11) allows us to write Ajcc, 0(2/2) m terms of the collinear factors at fixed values of 
their arguments, and thus to determine the 2/2 dependence of Ajcc, 0(2/2) for large negative y^\ 



AjCC,o(2/2) = AlCC, unren (M 2 ) exp 



In 



Cjcc, unren, W 

(0) + error, 



(18) 



with the error vanishing when y 2 —> —00. 

Exactly analogous steps show that the other unsubtracted quantities obey 



Bjcc,o(yi) = Bjcc, unren 

(M 2 ) exp 



2/i - In 



y/2k~ L 
M 



Cjcc, unren, W 

(0) + error, 



(19) 
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and 



^00,0(2/1,2/2) = CjCC, unren, 1 ty(0) exp 



(2/1 - 2/2 )#u 



Cjcc,w(0) + error. 



(20) 



The collinear factors associated with the Wilson lines are the same for all the Wilson lines. This is shown simply by 
exchanging the + and — coordinates. 

Inserting all of these results into AqBq/Sq shows that the dependence on y\ and 1/2 cancels, as do the factors of 
Cjcc,w- Therefore, the limit y\ —> +00 and y2 ^ —00 exists. It is just 



lim 



Ajcc,o(2/2)£jcc,o(2/i) 
Sjcc,o(z/1j2/2) 



AjCC, unren(M 2 ) B JC C, unren (^ 2 ) exp 
2 i unren (CO -^unren (Cb), 



m ^unren 



(21) 



given that 2k\k B = Q 2 and (a(b = Q 4 - The last result is exactly as in the CSS factorization formula, with the JCC 
definition of the collinear factors. 



B. Ao, Bo and So in the EIS method 



For the EIS method, we again apply the JCC factorization again to Ao, but with a changed regulator for the 
Wilson line. Instead of the non- light- like Wilson line of rapidity y<i that is used in (11), EIS use a ^-regulator. The 



steps for deriving factorization equations like (11) depend on the wide separation of rapidity between the hadron 
of momentum pa and the Wilson line in Aq. Thus the structure of factorization is insensitive to the method for 
regulating the Wilson line (so long as it is a method consistent with factorization), and the proof continues to apply. 
The factorization formula corresponding to (p~8|) is then 



^EIS,o(^ + ) = ^JCC, unren(2(/c+) 2 



-2yn 



Ceis, unren, W 

(y n — In 8 + ) + error, 



(22) 



where now the error vanishes as £ + — >• 0. In obtaining this equation, we have again used dimensional regularization. 
This regulates any soft and collinear divergences associated with massless quark and gluon lines, and therefore we 
have removed the A-regulator that EIS use. The collinear factor Ajqc is the same as in (11), because it is associated 
with the same collinear subgraphs associated with the hadron. But the different regulator for the Wilson line implies 
that the Cw factor is different. 

This factor, Ceis, unren, w(Vn ~ hic) + ), depends only on the combination y n — Yn.5 + but not on y n and 5 + separately, 
because of boost invariance. To see this result explicitly, we examine the relevant Wilson line propagators: 



(23) 



Suppose we simultaneously change y n to y n + Sy and S + to J + e 5y , so that y n — In 8 + is unchanged. Then we change 
variables on all loop momentum integrals in Ceis, unren, w by corresponding boost factors: k + to k + e 5y , and k~ to 



k~ 



-Sy 



This has effectively implemented a boost, and C is boost invariant, so its value is unchanged. 



In accordance with the JCC factorization method, Ceis 



: w is defined by 



Ceis, unren, w {yn ~ In 5+) = lim C E is,o (y[ ~ In 5+) 



^JCC,o(2/n,2/ 2 ) 



^JCC,o(2/i,2/ 2 ) ^JCC,o(2/i,2/n) 



C E is,o(+oo - ln(5 + ) 



^JCC,o(2/n, -OO) 



5'jcc,o(+oo,-oo) Sjcc,o(+oo,2/n) 



(24) 



The bare factor is defined with one Wilson line regulated by the EIS method and one regulated by the JCC method. 

The CS equations are the same as before, since they arise from the same 5jcc,o factors under the square root. 
Hence, instead of (18), we get 



-4EIS,o(£ + ) = AjcC, unren (M 2 ) exp 



m ^ -In — ] K u 

M M 



Ceis, unren, W 

(M) + error. 



(25) 
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This and the corresponding factorizations for Bq and So have errors that vanish as 5^ — >• 0. The collinear factor for 
the Wilson line changed from that in (11) and (18) because of the changed regulator. Note that the reference scale 
for the Ceis, unren,w(M) in Eq. (25) is now a mass M, since the regulator is a mass scale £ + . But the A factor is the 
same as in the JCC method. 

The formulas for Bq and Sq are 



#eis,o(£ ) = Ceis, unren,^ (- In 5 - y n ) B JC c, unren(2(/c B ) 2 e 2 ^) + error 



= Ceis, unren,w(M) exp 

SeIS,o(^ _ ,^ + ) = Ceis, unrei 

Ceis, unren,w(M) exp 



In 



V2k~r 



In 



Bjcc, 



M M 

^ w (- In 5~ - y n ) Ceis, unren.w (y n ~ In 5+) + error 
M 2 



,(M 2 ) 



In 



5-5+ 



Ceis, (M) + error, 



- error, 



(26) 



(27) 



Hence, the limit of AqBq/Sq is the same as in the JCC method: 



lim 



^EIS,o(^ + )^EIS,o(^~) 

S E is,o(^",^ + ) 



^-JCC, unren 

(M 2 ) Bjcc, unren (^ 2 ) GXp 



In , r n> I -^unrcn 



\/-') 



(28) 



As before, we have a cancellation of the collinear factors for the Wilson lines and of the dependence on the regulator 
parameters. Errors go to zero in the limit that the regulator is removed. Hence, we get the same finite result with 
the EIS method as with the JCC method. 

The above derivation used the JCC organization of factorization to express the combination AoBq/So of bare 
quantities in the EIS method for the Drell-Yan cross section in terms of A and B in the JCC scheme. But this same 
product is the product of parton densities in the EIS scheme. Hence from Eq. (28), it follows that the product AB is 



the same in both methods. The only question now is whether the individual factors A and B are the same. 

The only caveat in the above derivation arises because the derivation of factorization uses gauge invariance heavily, 
at least beyond simple one-loop graphs. But when Wilson lines in gauge- invariant operators are replaced by regulated 
Wilson lines in the EIS method, the resulting operators are not gauge invariant. This could potentially cause some 
problems. It needs to be investigated whether any problems actually occur. It is likely that any resulting corrections 
are suppressed by a power of the regulator parameters and will not affect the final results. 
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IV. THE JCC METHOD AND THE FORMULA A = A />/3o 
The unrenormalized collinear factors in the JCC method are defined by 



IjCC, unren (GO = Um A JC C,o(^2) J ^ £jCC,0 (^1 , 2/n) ^ 

y^+°° V ojcc,o(2/i,2/2) ^JCC,o(2/n,2/2j 



^JCC, unren (GO = Um SjCC.ufaO J c gjCC »o(^^ 2 ) ( 30 ) 

V 6 JCC,0(2/1j2/2) 6jCC,0(2/lj2/n) 

where ("a and are defined by Eqs. Q and (J5|. Evidently the product of Ajqc unren and 5jcc, unren is exactly (the 
appropriate limit of) AqBq/Sq, a result we have already seen in Eqs. (21) and (28). 

In these formulas, the limits y\ — » -f oo, 2/2 — oo can be taken independently. But the rate at which the rapidities 
go to infinity can also be coordinated. Let us do this so that the rapidity intervals are the same in both of the S 
factors that have a y n argument, i.e., let us set 

yi =y n + Ay, y 2 =y n - Ay. (31) 

Then we take the limit Ay — >> oo. Now, 

SjCC.oCz/l, Vn) = SjccAVn + A 2/, Vn) = ftcC,o(l/n, 2/n ~ = #JCC,o(2/n, 2/2) , (32) 

where the middle equality is obtained by exchanging the + and — coordinates, by applying charge conjugation to 
exchange the color charges of the Wilson lines, and by then applying a boost. After that, we apply UV renormalization. 
We therefore obtain 

Ajcc(Ca) = lim jcc, 0(^/2) x renorma 2i za ti on f ac tor, (33) 

11" -2 V #100,0(2/1,2/2) 

Bjcc(Cb) — li m — JCC, 0(2/1) x ren0 rmalization factor, (34) 
ll^tZ V^JCC,o(2/i,2/2) 

where we require the limits to be taken with ^(2/1+2/2) = Vn fixed. We can also write these with the limit of a single 
variable: 

Ajcc(Ca) — lim — JCC, 0(2/2) x renorma ij za ^i on f ac tor, (35) 

y^-°° V^JCC,o(22/n - 2/2,2/2) 

Bjcc(Cb) = lim Bjcc, 0(2/1) x renorma ii za ti on f ac tor. (36) 

y^+°° y/Sjcc,o(yi> 2 yn ~ 2/i) 

These formulas are simpler than the original JCC formulas, because they have two factors instead of four, and they 
are of the form of the EIS definitions except for the method of regulation of the light-like Wilson lines. 



V. EIS TMD PDFS HAVE A C PARAMETER AND ARE EQUAL TO THE JCC PDFS 

The EIS PDFs are defined by 

Aeis = lim — eis,o(^ ) x renormg Qj za tj on f ac tor, (37) 
^S EIS , (S-,5+) 

Beis — li m — ^eis,o(^ ) x jjy renorma ii za ^i on f ac tor. (38) 
V5 E is,o(^-^ + ) 

Compared with EIS's actual definition [3, Eq. (2.13)], we have changed the names of factors, we have made explicit 
the regulator parameters, and we have omitted the Fourier transform needed to get the TMD PDFs in transverse 
momentum space. 

In their statement of the definitions, EIS did not originally specify the relative rates at which £ + and 5~ go to 
zero. However, in their actual calculations — see [3, Sec. 3.1], they set J + = S~ . In fact, this is not a boost invariant 
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condition, and, as we now show, it suffices to keep the ratio 5 + /5 fixed as 5 + and 5 go to zero. From the perspective 
of factorization alone, any choice of fixed 8 + / 'S~ gives a distinct but equally valid definition in the limit of £ + , S~ — >• 0. 



The choice of 5 + = 5 , plus the choice to work in a frame where k~X 



Q, is what allows EIS to recover 



logarithms of Q 2 in their explicit calculation (e.g. [3j Eq. 3.10]), despite the apparent absence of the scale Q 2 inside 
their TMD PDF definitions. 



We apply (25) and (27) to (37), and apply UV renormalization, to get 



A E is = lim Ajcc(M 2 ) exp 



In 



V2kj 
M 



-In — 
2 5+ 



K 



(39) 



Evidently, the ratio S /5 + finite and nonzero for the limit to be well-defined. Equation (39) is exactly equal to the 
JCC PDF, Ajcc, with 



5+ 



(40) 



Then, in exact correspondence to the JCC scheme, we define 

U = 2{k + A fe-^=2{k\? 5 -, (41) 

and A E is depends on £a- Thus, the 5+,5~ regulation scheme of Ref. [3 is exactly equivalent to the JCC scheme if 
we identify the rapidity parameter as y n = In ^8+ /8~ . 

An exactly corresponding equation applies to the other PDF: 



^eis(Cb) 



lim Bjcc(M 2 ) exp 



In 



= £jcc(M 2 ) exp 



In 



M 



V2k} 
M 



Vn\K 



K 



with just an exchange of the roles of the + and — coordinates. Now we have 

C B = 2(fc^) 2 e 2 ^=2(fc^) 2 ^. 



(42) 



(43) 



The limit taken by EIS with 8 + = S~ simply corre- 
sponds to the case y n = 0, for which (a — (b — Q 2 , 
provided that one uses a frame where k\ — kg. Since 
the EIS PDFs have a £ argument, and since they equal 
the JCC PDFs, they must also obey a CS equation. In 
applications, the auxiliary parameters are to be evolved 
to values that give good or optimal accuracy for the per- 
turbative calculations that are relevant for a given pro- 
cess at a given hard scale. For the hard scattering, this 
entails setting \i and vT to Q (or to within a a finite 
factor of Q). In this sense the two auxiliary parameters 
can be treated as tied together. (In fact only the product 
CaCb of the two (s is relevant.) 

Although EIS [3 claim not to use the CS equation - 
see e.g., their Sec. 10 — they use an equivalent result, 
obtained in their Sec. 5 for the Q-dependence of their 
TMD pdfs. They label this result resummation, but the 
result is exactly equivalent to the use of the CS equation, 
the RG equation that is already part of their formalism, 
and the small-6 expansion of TMD pdfs. The results are 



the same as in the CS method (in the JCC version), as 
can be seen, for example, by comparing their Eq. (5.7) 
with Sec. 13.13.1 of Ref. [5]. EIS's restrict their results to 
the region A <C qr <C Q. However, with the exception of 
the small-6 expansion, the results also apply for all small 
qT as well, as is proved in Ref. [5]. 

Furthermore, EIS use here that the logarithm of the 
TMD PDFs are linear in ln(Q 2 //i 2 ). (See [3, Eq. (5.4)] 
and the discussion that follows.) In fact, this linearity is 
a trivial consequence of the CS equation. 

With the identifications in Eqs. ( 39p3 ), EIS find a one- 
loop anomalous dimension [3, Eq. (3721,3.22)] that is the 
same as the JCC one. 

EIS, in a recent preprint [7 , advocate setting y n = 
(£+ = S~) from the outset. While this is legitimate, 
we find it non-optimal for the following reasons. First, 
fixing a particular value for y n amounts to fixing a choice 
of reference frame with respect to boosts along the beam 
axis. In the case of y n = 0, it is the rest frame of the 
Drell-Yan lepton pair. This choice treats the hadrons 



11 



asymmetrically because, although it corresponds to the 
parton rest frame, k\ = kg = Q, it is not the rest frame 
of the incoming hadrons, so the fractional momenta xa 
and xb will in general be very different, so it does not 
correspond to a particularly natural frame with regard 
to the TMD pdfs. More importantly, setting y n to zero 
does not eliminate the CS evolution, but rather shifts it 
entirely into the dependence on k\ in Eq. (41) and kg 



in Eq. (43). In fact, CS evolution is boost-independent, 



which is a particular advantage of the JCC approach, so 
we view it as preferable to leave the frame unspecified in 
the implementation of evolution. 



VI. DIFFERENT VERSIONS OF MS 
RENORMALIZATION 

However a possible problem is that calculations of the 
hard part do not agree: Compare [19, Eq. (6)] with [3j 
Eq. (7.5)], where the tt 2 terms are different by a factor 
of 7/6. We suspect that this is due to a difference in 
the definitions of the MS scheme for UV renormalization. 
JCC [5 , Sec. 3.3.2] defines one- loop counterterms to have 
the structure 



S e x sum of poles in e, 



where 



r(i-, 



(44) 



(45) 



A more commonly used definition would replace S e by 
(47re _7s ) e , which differs by a term of order e 2 . Nor- 
mal UV renormalization only needs at most one factor of 
1/e per loop, and the change in the definition of the MS 
scheme is irrelevant for such quantities. But this is no 
longer true when there are UV divergences with a double 
pole 1/e 2 per loop, as in TMD pdfs with the EIS and JCC 
definitions. The modified definition Eq. (45) was made 



specifically with this situation in mind, since the factor 
S e is common to all one-loop transverse-momentum inte- 
grals in dimensional regularization. The use of (45) com- 



pletes the simplifications of the results of one-loop cal- 
culations provided by the MS scheme. We have checked 
that the difference in the coefficient of tt 2 is consistent 
with this change in the definition of the MS scheme. How- 
ever, not enough details are given in [3] for us to verify 
definitively that this is the correct explanation of the dif- 
ference between the reported hard-scattering coefficients. 



VII. SUMMARY 

In this article, we have shown that two seemingly dif- 
ferent methods of formulating TMD-factorization — one 
originating in ideas from SCET, and one based on sub- 
tractive methods of pQCD — are actually equivalent. 
That the same conclusions are found from two seemingly 



different approaches provides important support for the 
general validity of the TMD-factorization formalism. 

To clarify this equivalence, and to understand the 
origin of the superficial disparity between the two ap- 
proaches, it is necessary to comment on some differences 
between what we have referred to as the EIS method 



in Eq. (10) and the original presentation of the EIS for- 
malism in Ref. [3 : We have naturally assumed it to be 
implicit in Ref. [3] that the limit £ + , S~ — >• for the reg- 
ulators of the Wilson lines is to be taken, although this 
appears not to be explicitly stated. Now the calculations 
in Ref. [3] rely on separate regulators A + and A - for 
the quark lines, as well as regulators £ + and 5~ for the 
Wilson lines. The two kinds of regulator were tied by 
the conditions A + = S^Pg and A - = S~p\, so that in 
massless calculations, an appropriate soft approximation 
is valid. The A regulators cut off infrared divergences 
associated with the quark lines. If some other cutoff is 
imposed, for example, by nonzero masses, then the A 
regulator can be removed without being tied to the S 
regulator. In our general discussion, we have assumed 
that some other infrared cutoff is present throughout, 
and we have set A + = A - = from the outset. We note 
that the consistency of ^-regulator formalism is derived 
in Ref. [3] from an examination of one-loop graphs with 
on-shell external partons. It is not clear to us how com- 
plete the derivations Ref. [3] are beyond that order. For 
the purposes of this paper, we have taken as our starting 
point the assumption that overall factorization with the 
^-regulators remains consistent in a general derivation at 
arbitrary order. 

We made one other explicit extension to the EIS defi- 
nition. EIS [3 imposed the condition £ + = S~ . But this 
is boost dependent. We have instead allowed the ratio 
5 + /5~ to remain unspecified in Eq. (10) (its variation is 



responsible for the CS equation), and have observed that 
the EIS definition continues to work in this cas^fl 

Given the consistency assumption, we find not only 
that the methods (EIS and JCC) are equivalent, but 
that the different definitions give exactly equal TMD 
pdfs, with the exception of different definitions of the 
MS scheme, as explained around Eq. (45). A criticism 



that is sometimes made of the CSS formalism is that 
it is prohibitively difficult or complicated for practical 
use, particular in the context of probing non-pert urbative 
structure. However, we believe that the most recent ver- 
sion, presented in Ref. [5], has eliminated the main is- 
sues responsible for such difficulties, and now provides 
a maximally user-friendly TMD formalism. At least in 
the manual one-loop calculations we have performed, the 
length of the calculations is little different. The remain- 
ing complexities are those of ordinary Feynman graph 
calculations. 



4 We thank Markus Diehl for bringing to our attention the issues 
discussed in this and the previoius paragraph. 
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We conclude that since approaches to TMD factoriza- 
tion that are apparently widely different nevertheless give 
the same TMD pdfs, this supports the idea that there is a 
particularly natural kind of preferred definition of TMD 
pdfs in QCD. This gives good current and near- future 
prospects [20lf24] for applications of TMD-factorization 
in a broad variety of phenomenological studies. There 
remain differences in how the different groups analyze 
the non-perturbative large-6 region in the solution of the 
evolution equations. We leave the analysis of these dif- 
ferences to future work. The TMD pdfs themselves, once 
given an operator definition, are unambiguous entities in 
QCD. 

However, for precise phenomenological work there 
needs to be agreement on the definition of the MS scheme; 
differences in its definition matter for TMD pdfs even 



though they do not affect simpler quantities (where the 
UV divergences give only at most one factor of 1/e per 
loop). But this issue is not a matter of principle. 
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